Stringent Numerical Test of the Poisson Distribution for Finite Quantum Integrable 

Hamiltonians 
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Using a new class of exactly solvable models based on the pairing interaction, we show that it 
is possible to construct integrable Hamiltonians with a Wigner distribution of nearest neighbor 
level spacings. However, these Hamiltonians involve many-body interactions and the addition of a 
small integrable perturbation very quickly leads the system to a Poisson distribution. Besides this 
exceptional cases, we show that the accumulated distribution of an ensemble of random integrable 
two-body pairing hamiltonians is in perfect agreement with the Poisson limit. These numerical 
results for quantum integrable Hamiltonians provide a further empirical confirmation to the work 
of the Berry and Tabor in the semiclassical limit. 

PACS numbers: 02. 30. Ik (Sistcmas Integrables), 05.45.Mt (Quantum Chaos+Semiclassical chaos) 
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The concept of quantum chaos still lacks a clear defi- 
nition. The main ideas in this field have been obtained 
using the semiclassical approximation for quantum sys- 
tems havinga classical analogue. In a seminal paper, Bo- 
higas et al. [lj conjectured that the fluctuation properties 
of generic quantum systems, which in the classical limit 
are fully chaotic, coincide with those of random matrix 
theory (RMT). This conjecture is strongly supported by 
experimental data, numerical calculations, and analyti- 
cal work based on semiclassical arguments. For a generic 
quantum integrable system, Berry and Tabor Q showed 
that under very general conditions the spectral fluctua- 
tions in the semiclassical limit are well described by the 
Poisson statistics, i. e. the successive energy levels are 
uncorrelated. 

The analysis of spectral fluctuations provides an es- 
sential tool in the study of quantum chaos. Moreover, 
as RMT was introduced to explain the fluctuation prop- 
erties of many-body systems like the atomic nucleus, it 
is usually considered that the level statistics establishes 
a link between many-body systems and the semiclassical 
picture. Actually, if the semiclassical limit is not valid, 
the comparison of the system spectral fluctuations with 
those predicted by RMT is the main criterion to decide 
whether the system dynamic is regular or chaotic. When 
spectral fluctuations fall very near the RMT predictions 
the quantum system is considered fully chaotic. On the 
contrary, if they follow closely enough the Poisson statis- 
tic, the system is considered regular. Therefore, the con- 
cepts of chaotic and regular quantum systems are not 
well established since they rely on results that have only 
been checked in the semiclassical limit. 

The concept of integrability in Classical Mechanics is 
well defined after the work of Liouville in the 19th cen- 
tury. A classical Hamiltonian system is integrable if it 
has a complete set of independent integrals of motion 
commuting with respect to the Poisson brackets. The 
total number of integrals of motion should be half of the 
dimension of the phase space to assure completeness. In 



Quantum Mechanics, the concept of integrability is usu- 
ally derived from an extension of the Liouville's defini- 
tion. A quantum system is said to be integrable if it is 
possible to define a complete set of hermitian operators, 
the integrals of motion, that commute among each other. 
However, this definition has some deficiencies and ambi- 
guities. Various attempts to clarify the characteristics of 
these integrals of motions, mainly their functional inde- 
pendency, have not produced irrefutable answers (see, for 
instance, 0- [J 13 and references therein). In this work 
we shall use the following definition: a quantum system is 
said to be integrable if a set of as many commuting her- 
mitian operators as quantum degrees of freedom can be 
explicitly given, and the hamiltonian can be expressed as 
a function of these operators Q . This criterion, relying in 
the algebraic structure of Quantum Mechanics, does not 
directly refer to Classical Mechanics. Since this defini- 
tion requires that a "complete" set of quantum integrals 
of motion is explicitly given, it assures the existence of a 
common basis of eigenstates labelled by the eigenvalues 
of the integrals of motion. Moreover, the system is Ex- 
actly Solvable if the complete set of eigenstates can be 
found by algebraic methods. In this sense, the previous 
definition of quantum integrability is closely related to 
exact solvability. We will use this criterion to test some 
accepted properties of quantum integrable systems, like 
the Berry and Tabor conjecture, without taking into ac- 
count the classical limit. 

The level statistics of quantum integrable models has 
been analyzed in condensed-matter physics as well as in 
nuclear physics. Poilblanc et al. Q studied the spectral 
fluctuation by finding the energy spectrum for several 
one-dimensional finite lattice models, like the Heiscnbcrg 
model, the t — J model and the Hubbard model. Al- 
hassid and Novoselsky Q studied the quantal fluctua- 
tions displayed by the energy levels in the Interacting 
Boson Model of nuclei. Recently, the realization of a 
Poisson distribution has been suggested as a detector of 
new integrable quantum chiral Potts models Q . In each 
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case, the intcgrablc Hamiltonian is parameter indepen- 
dent, or it depends on a single free parameter like in the 
Hubbard model (the on-site repulsion U). The study of 
the spectral fluctuations has been carried out diagonaliz- 
ing a definite Hamiltonian in the largest possible Hilbcrt 
apace. In all cases of integrable quantum models, it was 
verified that the histogram of the near-neighbors level 
spacing could be well fitted by a Poisson distribution. 
These results provide a numerical support for the Berry 
and Tabor semiclassical demonstration H, Q, , but the 
quality of the statistics is poor due to the small num- 
ber of energy levels entering in the histogram. This is 
in contrast with the numerical testing of the Bohigas [| 
conjecture, assigning a Wigner distribution to any non- 
integrable quantum hamiltonian, in which large ensem- 
bles of two-body random Hamiltonians were considered 
(for a review see @). More recently, Bcnct et al. [Toj 
have studied an ensemble of integrable bosonic hamilto- 
nians whose members display GOE or GUE-like spectra 
with probability one. This anomalous behavior can be 
explained in terms of semiclassical mechanics. Despite 
of the fact that the systems under consideration were in- 
tegrable (in the semiclassical sense), the periodic orbits 
that fulfill the Einstcin-Brillouin-Kellcr (EBK) quatiza- 
tion condition explore huge regions of the phase space, i. 
e. they mimic a typical chaotic motion; thus, it is reason- 
able to find random matrix spectral fluctuations. There- 
fore, the whole ensemble can be considered an anomalous 
exception of Berry and Tabor conjecture. 

Trying to get more insight on these ideas, we study 
the level statistic of a new class of quantum intcgrablc 
models: the Richardson-Gaudin models; they are based 
on the pairing interaction |lll | and have a large number 
of free parameters which can be selected randomly. In 
particular, we will study the rational model for which 
the quantum invariants have the form 




+K°K°}, (1) 

where i labels the fl levels of a single particle basis, rji 
are fl free real parameters and g is the paring strength. 
The operators K + , K~ and K° are the SU(2) generators 
of the pair algebra in level i 

X? = i(oto i + o} 0f -l) , Kt= l -a\al t = {K7)\ 

m 

The operator Kf creates a pair of particles in time 
reversal states in the double degenerated level i. The 
three generator close the commutation algebra SU (2) 

[Kf, K+] = 5 U >K+ , [K+, K-] = 26 u ,Kf. (3) 

Based on this commutation algebra, it is straightfor- 
ward to check that the operators commute with one 



another for arbitrary values of g and the set of Q param- 
eters rj. The rational model, as well as other models, is 
fully integrable and exactly solvable (for the exact eigen- 
states of these models see [Tll|'). 

Once the free parameters inside the Ri operators are 
fixed, it is possible to find their complete set of common 
eigenvalues rf and eigenvectors \ip a >. Any function of 
the Ri operators defines a valid integrable Hamiltonian. 
In particular, linear combinations of the Ri operators 
produce quite general pairing hamiltonians 

i i 

9 E ^ {\ i K t*i + KTKj) . (4) 

Contrary to most of the intcgrablc Hamiltonians cited 
in the introduction, these new class of integrable Hamil- 
tonians depend on (2f2 + 1) free, real and independent 
parameters; combinations of higher rank quantum invari- 
ants would give larger sets of free parameters. Within 
this very large parameter space it is worth to explore the 
possibility of finding chaotic spectra, contradicting the 
Berry and Tabor hypothesis. Therefore the usual dis- 
tinction between chaos and regularity by means spectral 
statistics would become doubtful. 

To shed some light on this question we have tried to fit 
several chaotic spectra using the class of Hamiltonians Q 
for different values of Q and different particle pair num- 
bers N. These chaotic spectra were obtained by diagonal- 
izing a random matrix with the appropriate dimension. 
The dimension of the Hilbcrt space for a system with £1 

fn\ 

levels and N fermion pairs is given by D = I J . Thus, 

is usually much larger than the number of Hamiltonian 
parameters 20 + 1; therefore it is impossible to obtain 
an exact replica of the actual random matrix spectrum, 
and we can only get an approximation to this spectrum. 
In all the cases the best parameter set leads to a Hamil- 
tonian whose spectrum shows Poisson level fluctuations, 
as predicted by Berry and Tabor. 

Nevertheless, as we have commented above, it is still 
possible to consider more general Hamiltonians using the 
Ri operators. These Hamiltonians involve many body 
forces represented by combinations of higher rank R op- 
erators. Knowing the dimension D of the Hilbcrt space 
for a system with Q levels and N fermion pairs, we pro- 
pose the following class of intcgrablc Hamiltonians with 
many body forces. 

il>i2>--->iN 

If we fix the rji parameters and the value of g inside the 
Ri operators entering in J5J , the number of free param- 
eters Si 1 _i 2 _...j N in the hamiltonian is precisely equal to 
the dimension of the Hilbcrt space D. Therefore, the d 
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eigenvalues of a non-integrablc pairing Hamiltonian with 
a typical chaotic spectrum can be exactly fitted with the 
Hamiltonian ijjjj by solving a linear set of equations for 
the D unknowns £ii,i 2 ,---,ijv 

We have considered several cases with (fi, N) = (10, 6) 
(11, 6), (12, 6) and (13,6); the corresponding dimensions 
of the Hamiltonian matrices are D = 210, 462, 924, 1716 
respectively. In each case we have been able to fit a 
chaotic spectra with the appropriate dimensionality ob- 
taining an exception to the Berry and Tabor result. In 
other words, we have been able to obtain spectral fluctu- 
ations of Wigncr-Dyson type at all the energy scales for 
quantum integrable hamiltonians. Actually, exceptions 
to this rule were already known: for example, Crehan 
|l2j proved that any spectral sequence obeying a certain 
growth restriction is the quantum spectrum of an equiv- 
alence class of classically integrable nonlinear oscillators. 
Our result, however, is quite more enlightening because 
it provides an example of a full and realistic quantum 
system where the semiclassical result of Berry and Tabor 
does not apply. 

Our findings rise the question of how stable are these 
solutions against small perturbations of the Hamiltonian 
parameters within the parameter space and whether fi- 
nite size effects may affect these properties. To study 
these two points the Hamiltonians obtained in the pre- 
vious fit are perturbed as follows. Each parameter 
E *i,ia,— ,»jv m © i s replaced by 

e — >e =e(l + <j>\), (6) 

where A stands for the perturbation strength and <^> is a 
phase chosen at random. Notice that the new Hamilto- 
nian is also a combination of the Ri operators and there- 
fore it is integrable. 

Before we proceed to analyze the spectral fluctuations 
of the perturbed Hamiltonians, their spectra must be un- 
folded. For any quantum system the level density p(E) 
can be separated into a smooth part p(E) and a fluc- 
tuating part p(E). The former gives the main trend of 
the level density and the later characterizes the spectral 
fluctuations. Similarly, the cumulative level density, that 
gives the true number of levels up to energy E, 

N(E) = f dE'p(E'), (7) 

can be separated into a smooth part p(E) and a fluctu- 
ating p(E) part, i.e., 

N(E) = N(E)+N{E). (8) 

Actually, level fluctuation amplitude s are modulated by 
the value of the mean level density p{E); therefore, to 
compare the fluctuations of different systems, or even the 
fluctuations of different parts of the same spectrum, the 
level density smooth behavior must be removed. This is 



1 




' -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 

LoglOW 

FIG. 1: Variation of the Brody parameter ui as a function of 
the perturbation strength A for (fi, N) = (10, 6) (11, 6), (12, 6) 
and (13, 6). 



done by means of a transformation, called unfolding [l3j , 
which consists in mapping the energy levels Ei onto new 
dimensionless (i = N(E). Then, the nearest neighbor 
spacing sequence is defined by 

Si = ( i+ i-(i, i = l,-- ,N-1. (9) 

For the unfolded levels p(e) — 1 and (s) = 1. In general 
this is a difficult task for systems where an analytical 
expression for the mean level density is not known. This 
is the case of the hamiltonian ensembles introduced in 
the present Lette r. The refore, in order to obtain a good 
approximation to N(E) we have performed a least square 
fit to Chebyshev polynomials. 

The spectral properties of the perturbed Hamiltonians 
can be analyzed by means of different statistics. The 
most simple is the nearest neighbor spacing distribution 
P(s). It allows to quantify the chaoticity of the sys- 
tem in terms of a single parameter oj by fitting the P(s) 
histogram to the Brody distribution P(s,uj) [14j using a 
least-squares method. The resulting distribution interpo- 
lates between the Poisson limit (uj = 0) and the Wigner 
limit (ui = 1). It would be also possible to analyze the be- 
havior of the eigenstates as a function of the perturbation 
strength. This could be done calculating the information 
entropy Ih or the localization length Ih |15|. However, 
we can advance that the Hamiltonian eigenstates will re- 
main unaltered under the influence of the perturbation 
(JHJ) which is defined on the coefficients of the combina- 
tion of integrals of motion but does not modify them. 
Actually, for any integrable system the eigenstates are 
completely defined by the quantum integrals of motions, 
and the Hamiltonian can be expressed as a function of 
these operators. Therefore, all Hamiltonians obtained us- 
ing different functions of a given complete set of integrals 
of motion will have the same eigenstates. 

Fig. n displays in a semi-logarithmic scale the Brody 
parameter u> as a function of the perturbation strength A 
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for the four (f2, TV) examples we have considered before. 
In all the cases a very small perturbation is enough to 
drive the system to the Poisson limit. For the smallest 
system we considered (D = 210) a perturbation A ~ 
10~ 5 is enough to obtain poissonian spectral fluctuations, 
while for the largest system (D = 1716) three order of 
magnitude less are required. Clearly, the trend is that 
larger systems require smaller perturbations. 

On the light of these results we conclude that it is nec- 
essary to consider integrablc hamiltonians with many- 
body interactions in order to obtain a chaotic-like energy 
spectrum. However, small perturbations within the in- 
tegrable space of parameters restore the poissonian like 
spectrum. We conjecture that for very large f2 and N val- 
ues a chaotic spectra would correspond to isolated points 
in the parameter space and that a infinitesimal perturba- 
tion within this space would immediately drive the sys- 
tem to a Poisson distribution. 

In some cases it is possible to introduce a suitable ran- 
dom matrix model that describes the behavior of spectral 
fluctuations as the system evolves through the parameter 
space. This has already been done for the metal-insulator 
transition or for the order-disorder transition in Quantum 
Hall systems |l6j . However, in the prersent case, this ap- 
proach seems to be more complicated since the "Wigner 
phase" corresponds to isolated points in the parameter 
space. 

In view of the previous results, one would expect that 
a more physically sound family of pairing Hamiltonians 
with two body forces like l@J should give rise to a clean 
Poisson level statistics. In order to verify the correct- 
ness of this statement we have studied an ensemble of 
Hamiltonians l@J that we shall call "Two Body Pair- 
ing Random Ensemble" (TBPRE). We selected the case 
of (fl, N) = (13,6), corresponding to largest dimension 
D = 1716. For the sake of simplicity, the rji param- 
eters and the strength g were given fix values and the 
Si coefficients were selected randomly. The quantum in- 
variants Ri arc independent of the Ei parameters and 
thus, their corresponding eigenvalues stay the same for 
the whole ensemble. In our calculations we have used 200 
ensemble members. 

The short and long range spectral correlations of this 
ensemble have been analyzed by means of the usual level 
statistics distribution P(s) and by the rigidity A 3 (L) re- 
spectively. In the Poisson limit, characteristic of a reg- 
ular system, the nearest neighbor spacing distribution 
behaves as p^«s«°»( s ) = exp (_ s ) an d Af OISSOIl (L) = 
i/15. Fig. |3 compares the nearest neighbor distribution 
P(s), calculated numerically for our ensemble with the 
expected Poisson limit. We present the results in normal 
as well as in semi-logarithmic scales to enlarge the small 
spacing region, which has been shown to present some 
deviations from the Poisson limit [13 . ^ can be seen 
that the histogram and the theoretical curve match per- 
fectly. Fig. |3 shows the calculated A 3 (L) and compares 
it to the Poisson limit L/15. The agreement is almost 
perfect up to L ~ 40; for larger L values the A3 shows 




FIG. 2: Nearest neighbor spacing distribution, P(s), for 200 
TBPRE members. The dashed curve correspond to the Pois- 
son limit. 



an slight upbending from the Poisson straight line. The 
calculation of the rigidity is quite sensitive to unfolding 
procedure. When the mean level density is not accurate 
enough known, the unfolding method will introduce ac- 
cumulated errors that eventually give rise to a spurious 
increase of the A 3 (L) for large L values 01- The results 
shown in these two figures make it possible to conclude 
that the TBPRE spectral fluctuations are very well de- 
scribed by the Poisson statistic. Actually, Figures 2 and 
3 constitute the most precise numerical verification of the 
Berry and Tabor theoretical proof due to the fact that we 
were able to accumulate statistics by using an ensemble 
of random integrablc Hamiltonians (TBPRE) which, to 
the best of our knowledge, it wouldn't be possible for any 
other intcgrable model. 
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FIG. 3: A 3 (L) statistic calculated for 200 TBPRE members. 
The dashed line represents the Poisson values. 



To support in a more precise way this conclusion we 
consider again the p(s) statistics. As this statistics is less 



5 



sensitive to the unfolding procedure than the A3, small 
deviations from the theoretical Poisson prediction can be 
connected to the actual characteristics of the system dy- 
namics. To avoid any effects related to the bin size in 
the P(s) histogram, we will use the accumulated nearest 
neighbor spacing distribution 7(s) = J Q S P(s')ds' , which 
in Poisson limit it is given by j Pmsson (g) = l — cxp (— s). 
We define a "distance" between the calculated n(s) dis- 
tribution and the Poisson limit as 



S 2 = \n(s) -n Pmsson (s)\ ds. (10) 
Jo 



We consider four TBPRE with {Sl,N) = (11,6), 
(12,6), (13,6) and (14,6). The matrix dimensions for 
these ensembles are D = 462, 924, 1716 and 3003 respec- 
tively. In order to have approximately 2 x 10 5 spacings 
in each of the four ensembles, different number of mem- 
bers were chosen for each one of them. Fig. 0] shows the 
logarithm of the average distance S as a function of the 
matrix dimensionality D. The most relevant observed 
feature is that S decreases by an order of magnitude as 
the dimensionality increases from D = 462 to D = 3003. 
Moreover, the smooth and monotonous decrease of this 
function suggests that it goes to zero in the large D limit. 
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FIG. 4: Logarithmic plot of the average quadratic distance, 
5 2 , between the accumulated nearest spacing distribution ob- 
tained from the 200 TBPRE members and the Poisson limit, 
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given by n [s) = 1 — exp {—s). 



Summarizing, the use of a new family of fully intc- 
grable and exactly solvable pairing models with a large 
number of free parameter which can be selected at ran- 
dom allowed us to perform several stringent tests of the 
Berry and Tabor semiclassical proof. Based on the nu- 
merical results obtained, we conclude that quantum in- 
tegrable systems indeed follow a Poisson distribution of 
nearest neighbor level spacings for large enough systems. 
Exceptions to this rule can be found, but we showed that 
they are isolated solutions of high rank Hamiltonians and 
that they quickly decay to a Poisson distribution with an 
infinitesimal integrable perturbation. 



We thank O. Bohigas, P. Leboeuf and G. Sierra for use- 
ful discussions. This work has been supported by grants 
BFM2003-05316-C02-02 and BFM2000-0600 



[1] O. Bohigas, M. J. Giannoni and C. Schmit, Phys. Rev. 

Lett. 52, 1 (1984). 
[2] M. V. Berry and M. Tabor, Proc. R. Soc. London A 356, 

375 (1977). 
[3] St. Weigert, Physica D 56, 107 (1992). 
[4] S. Weigert and G. Muller, Chaos, Solitons and Fractals 

5, 1419 (1995). 

[5] W. M. Zhang and D. H. Feng, Phys. Rep. 252, 1 (1995). 
[6] D. Poilblanc, T. Ziman, J. Bellissard, F. Mila and G. 

Montambaux, Europhys. Lett. 22, 537 (1993). 
[7] Y. Alhassid and A. Novoselsky, Phys. Rev. C45, 1677 

(1992). 

[8] J. Ch. Angles d'Aur iac, J. M. Maillard and C. M. Viallet, 

|co^ -mat/020 5l"0l"1 
[9] V. K. B Kota, Phys. Rep. 347, 223 (2001). 
[10] L. Benet, F. Leyvraz and T. H. Seligman, Phys. Rev. E 



68, 045201(R), (2003). 
[11] J. Dukelsky, C. Esebbag and P. Schuck, Phys. Rev. Lett. 

87, 066403 (2001). 
[12] P. Crehan, J. Phys. A 28, 6389 (1995). 
[13] F. Haake, Quantum Signatures of Chaos, (Springer- 

Verlag, Berlin, 2001). 
[14] T. A. Brody, J. Flores, J. B. French, P. A. Mello, A. 

Pandey and S. S. M. Wong, Rev. Mod. Phys. 53, 385 

(1981). 

[15] F. M. Izrailev, Phys. Rep 196, 299 (1990). 

[16] P. Shukla, Phys. Rev. E 62, 2098 (2000). 

[17] G. Casati, B. V. Chirikov and I. Guarneri, Phys. Rev. 

Lett. 54, 1350 (1985). 
[18] J. M. G. Gomez, R. A. Molina, A. Relafio and J. Reta- 

mosa, Phys. Rev. E 66, 036209, (2002). 



